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Abstract. We consider an ensemble of classical particles coupled to a Klein- 
Gordon field. For the resulting nonlinear system of partial differential equa- 
tions, which we call the relativistic Vlasov-Klein-Gordon system, we prove the 
existence of global weak solutions for initial data satisfying a size restriction. 
The latter becomes necessary since the energy of the system is indefinite, and 
only for restricted data a-priori bounds on the solutions can be derived from 
conservation of energy. 

1. Introduction 

When considering the interaction of classical particles with classical or quantum 
fields various different situations arise: On the one hand one can consider the cou- 
pling of a single classical particle to a field. How this should properly be done for 
the case of a Maxwell field is a classical problem, cf. and the effective dynam- 
ics and asymptotics of such systems is an active field of research, cf. [0, ^, |l^, |f) 
and the references there. On the other hand, in kinetic theory ensembles of classi- 
cal particles are considered which interact by fields which they create collectively. 
There is an extensive literature on such systems, with particles interacting by non- 
relativistic, gravitational or electrostatic fields — the Vlasov-Poisson system — , by 
electrodynamic fields — the Vlasov-Maxwell system — , or by general relativistic 
gravity — the Vlasov-Einstein system. In all these systems the only interaction 
of the particles is via the fields which they create collectively, a situation which is 
sometimes referred to as the mean field limit of a many-particle system. 

In the present paper we consider an ensemble of particles which can move at 
relativistic speeds and interact by a quantum mechanical Klein-Gordon field. Let 
f = f{t,x,v)>0 denote the density of the particles in phase space, p = p{t,x) their 
density in space, and u = u{t,x) a scalar Klein-Gordon field; tgR, xgR'^, and 
€ denote time, position, and momentum respectively. The system then reads 
as follows: 

dtf + v-dJ-d^u-d^f^O, (1.1) 
d^u-Au + u = -p, (1.2) 

pit,x)= / fit,x,v)dv. (1.3) 
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Here we have set all physical constants as well as the rest mass of the particles to 
unity, and 

, (1.4) 

denotes the relativistic velocity of a particle with momentum v. We call this system 
the relativistic Vlasov-Klein-Gordon system. 

To our knowledge it has not yet been considered in the literature. Our motivation 
for initiating a study of this system is the following. In p a single classical particle 



coupled to a Klein-Gordon field is considered, and the system (1.1), (1.2), ( p_.3| ) is 
meant as a natural generalization of this to the many-particle situation. On the 
other hand, the system falls within the general class of nonlinear PDE systems like 
the Vlasov-Poisson or Vlasov-Maxwell system, and by studying it one may hope to 
learn more about the general properties of this important class of problems from 
mathematical physics. A major issue in this area, which we also focus on in the 
present paper, is the question of global existence of solutions to the corresponding 
initial value problem. 

The coupling in the system above is set up in such a way that the system is 
conservative; 

J y^l + \v\^fdxdv+^J [|9(Up + |a^u|2 + |u|2]da; + J pudx=:EK + EF + Ec 

(1-5) 

is conserved along sufficiently regular solutions. From experience with the related 
systems from kinetic theory mentioned above one knows that for global existence 
questions it is essential to derive a-priori bounds on the solutions from conserva- 
tion of energy. In this context the relativistic Vlasov-Klein-Gordon system poses 
the following specific difficulty: The interaction term Ec in the energy need not 
be positive. Therefore, one has to try to estimate this term in terms of the pos- 
itive quantities Ek and Ep- By Holder's inequality, interpolation, and Sobolev's 
inequality. 



pudx 



< 



\\pme/5\\ume<CEK{tndMt)\[ 



the details of a more general version of this estimate can be found below, cf. 



Lemma 4.1. The problem is that, after applying the Cauchy inequality, the right 
hand side in this estimate is of the same order of magnitude as the positive terms 
in the energy, and no a-priori bound on the solution seems to follow from conserva- 
tion of energy. This situation is similar to the gravitational case of the relativistic 
Vlasov-Poisson system where it is known that spherically symmetric solutions with 
negative energy blow up in finite time B . Our way out of this difficulty is to observe 

o o o 

that in the estimate above the constant C depends on ||/||i and ||/||oo, / being the 
initial datum for /, and if we require that it is smaller than an appropriate threshold 
then a-priori bounds on the solutions can be derived from conservation of energy. 
These bounds are such that one can pass to the limit along a sequence of global 
solutions to an appropriately regularized system, and this limit is a weak solution 
to the relativistic Vlasov-Klein-Gordon system. The details of these arguments are 
then similar to the corresponding ones for the Vlasov-Maxwell system cf. also 

HI. 

The paper proceeds as follows: In the next section we collect for easier reference 
some results on the linear, inhomogeneous Klein-Gordon equation. In Section 3 we 
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prove a global existence and uniqueness result for smooth solutions to an appro- 
priately regularized version of the Vlasov-Klein-Gordon system. We then proceed 
to derive uniform a-priori bounds on a sequence of such regularized solutions from 
conservation of energy for restricted data, and prove the existence of a global, weak 
solution to the original system. This is done in Section 4, where we also discuss 
some of the properties of these weak solutions. 

We conclude this introduction with some remarks on related results in the lit- 
erature; we restrict ourselves to results on the initial value problem in the three 
dimensional case. For the Vlasov-Poisson system global classical solutions for gen- 
eral initial data have been established in ||l6|, For the relativistic Vlasov- 
Maxwell system such global classical solutions are so far only known for small, 
nearly neutral, or nearly spherically symmetric data P, ^, |l^ . For general data 
the existence of global weak solutions was obtained in p| . When passing from the 
Vlasov-Poisson to the Vlasov-Maxwell system a lot of the difficulties with classical 
solutions of course arise from the different field equation which becomes hyperbolic 
instead of elliptic. However, it must be emphasized that already for the relativistic 
Vlasov-Poisson system where the only difference to the Vlasov-Poisson system is 
that a t; in the Vlasov equation is replaced by a u no global existence result is known 
for classical solutions with general initial data, and for the gravitational case where 
the energy is indefinite classical solutions can blow up in finite time as mentioned 
above, cf. |3|. An investigation of the fully relativistic gravitational situation, i.e., 
of the Vlasov-Einstein system, was initiated in |]l8) . 



2. The linear Klein-Gordon equation 

Although our notation is mostly standard or self-explaining we explicitly mention 
the following conventions: For a function h — h{t,x,v) or h = h{t,x) we denote for 
given t by h(t) the corresponding function of the remaining variables. By || .||p we 
denote the usual L^'-norm for p€ [l,c»]. The index c in function spaces refers to 
compactly supported functions. Throughout the paper the convolution denoted by 
* refers to the spatial variables. 

For easier reference we now collect some results on the linear, inhomogeneous 
Klein-Gordon equation 

dtU-Au + u = g (2.1) 
with a prescribed right hand side g^g{t,x). For t>0 and xgM^ let 
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47r t 4tt ' " y/t^-\x\ 

where S is the (5-distribution, H is the Heaviside function and Ji is the Bessel 
function of the first kind. 
Given initial data 

u{0)=uu dtu{0) = U2 (2.2) 



the solution of the initial value problem (2J), ( p.2[ ) can be written as 

u{t,x) —Uhom{t,x)+Uinhit,x), t>0, XGM^. 
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Here 



-i- ( U2{x-y)dSy-^ [ Mx-y)^^^dy 



with the abbreviation 



is the solution of the homogeneous Klein-Gordon equation ( p.lD with initial data 
(H), and 

Uinhit,x) := J J R{t- s,x~y)g{s,y)dyds 

7- / / 9{s,y)dSy-^-^ I I g{s,y)'^^^dyds 

47rJo J\x-y\=t-s t-S AttJq J\x-y\<t-s ? 

with the abbreviation 

i:^^{t-sf-\x-y\^ 

is the solution of the inhomogeneous Klein-Gordon equation with vanishing initial 
data. These formulas can be found in and They can be established by 

observing that the substitution w{t,xX) '■=u{t,x)eyip{—i(^) transforms (2.1) into a 
wave equation for w which can be solved in the usual way. 

If a G Nq denotes an arbitrary multi-index and d" the corresponding spatial 
derivative the above formulas imply the following estimate for the solution of (2.1), 
(2.2), which is not very sophisticated but good enough for our purpose: 

||5"u(t)lloo<C(l + <)4(||a"ui||oe + ||Va"wi||oo + ||5"u2||oo + ||9"g(t)||oo). (2.3) 



3. Global classical solutions of the regularized system 

The aim of the present section is to establish the following global existence and 
uniqueness result for a suitably regularized Vlasov-Klein-Gordon system: 

Theorem 3.1. Letd^C^{MP). Consider the regularized, relativistic Vlasov-Klein- 
Gordon system where the right hand side —p in is replaced by ~p*S. For 

initial data 0</gC^(R^), ui gC^{M.^) and U2 £C^{M.^) there exists a unique solu- 
tion {f,u) to the regularized system with f €C^{[0,oo[xM.^) anc? w G C^([0,oo[xR'^), 

o 

satisfying the initial conditions f(0)^f, u(0)—ui, dtu{0) —U2- 

Proof. We begin by recursively defining iterates /„ : [0,oo[xR^ ^ [0,oo[, u„ : 

o 

[0,oo[xM'^^R: Let fo{t,z) :— f{z), z—{x,v), and suppose /„ has already been 
defined. Let pn{t,x) :— J fn{t,x,v)dv and let u„ denote the solution to the Klein- 
Gordon equation ( |2.l| ) with right hand side —pn*S and initial data ui, U2, cf. 
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the previous section. Denote by Zn{s,t,z) — {Xn,Vn){s,t,x,v) the solution of the 
characteristic system 

X — V, 

V = -dxUn{s,x) 

with initial datum Zn{t,t,z) = z. The next iterate is then defined as fn+i{t,z) := 

o I 1 

f{Zn{Q,t,z))^ i.e., as the solution of the Vlasov equation (O) with u replaced by 

o 

Un and initial datum /. 

The flow of the characteristic system is measure preserving, hence = 
||p„(t)||i = ||/||i. This implies that for any aSN^, ||5"(/9„(i)*<5)||oo is bounded, 
uniformly in t > and n € N. Let T > be arbitrary. In what follows constants 
denoted by C may change from line to line and depend on the initial data, on the 
regularization kernel 5, and on T, but never on nSN. By ( |2.3| ), ||9"M„(t)||tx3 < C' 
for \a\ <2 and t£ [0,T]. Hence there exist R>0 and P>0 such that fn{t,x,v) = 
if >P or \x\>R and t€ [0,T], in particular, /„(t) eC^(R'^). Using these bounds 
it is now easy to show that (/„)„gN is a uniform Cauchy sequence on [0,T] xM^: 
By definition, 

||/n+lW-/n(i)||oo<C||Z„(0,i,.)-Z„_i(0,t,.)||oo. 

Using the characteristic system and the bound on 9^u„(t) we obtain, abbreviating 
Zn{s,t,x,v) by Zn{s), 

\Xn{s)~XnMs)\< j \Vn{T)-Vn-l{T)\dT 
J s 

\Vnis)^Vn-i{s)\ <cj^ |X„(T)-X„_i(r)|dr 

+ / \\dxUn{T) - dxUn-l{T)\\ocdT, 
J s 

hence 

|Z„(0)-Z„_i(0)| <c / I {T)\dT+ / WdxUnir)- dxUn-l{T)\\oodT. 



By Gronwall's inequality, 

\\fn+l{t)- fn{t)\\oc<C / \\dxUn{T)-dxUn-l{T)\\oodT. 

Jo 

On the other hand, the formulas for u„ from the previous section, the fact that we 
have regularized the right hand side of (1^), and the above bounds on the support 
of /„ imply that 

WdxUnir) - dxU,i^i{T)\\oc < C||p„(t)-p„_i(t)||i 

< C* / / \fn{T,X,v)- fn-l{T,X,v)\dvdx 



< C\\fn{T)-fn-lir)U. 

Combining these estimates we obtain 

1 1 /«+ 1 (i ) - /« (i ) 1 1 oo < C / V« (t) - /„- 1 (r) I loo dr. 
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Thus, by induction, 

||/n+lW-/nWI|oo<C— - 



SO (/n)nGN IS Uniformly Cauchy on [0,r] xR®. The same is true for {pn)neN and 
(9"w„)„gN for \a\ <2. That the uniform Umit of the iterative sequence has 

the required regularity and is the unique solution of our initial value problem on 
[0,r] follows, and since r>0 was arbitrary the proof is complete. □ 
In the next section we want to obtain a global weak solution of the relativistic 
Vlasov-Klein-Gordon system as a limit of a sequence of solutions to systems regu- 
larized with 5n 's that converge to the ^-distribution. To do so we will need energy 
bounds on these regularized solutions, but since we have modified the system en- 
ergy conservation takes a somewhat different form from what was stated in the 
introduction: 



Lemma 3.2. Let deC^{W) be even, 5 = d*d, and ui, U2€CciW). Let 0</G 
C^(M^) and let {f,u) he the unique solution to the regularized system according 

o 

to Theorem \3.\ with initial conditions f{0)~f, u(0)—ui*S, dtu(0) — U2*S. Let 
p{t,x) — J f{t,x,v)dv and denote by u the unique solution to the initial value prob- 
lem 

dfU — Au + u~—p*d, 
u{0)—ui*d, dtu{0)=U2*d. 

Then 

E:= j ^/l + \v\^fdxdv+]^j [\dtu\^ + \d^u\^ + \u\^]dx + j pudx =: Ek + Ep + Ec 
is constant in t. 

Proof. Using the Vlasov equation and integration by parts we obtain 

— = - / dtpudx. 

Also, 

— = — y dtup^ ddx ~ — J dt{u*d) pdx] 

for the last equality observe that d is assumed to be even. The convolution u*d 
satisfies the Klein-Gordon equation with right hand side —p*S and initial conditions 
{u*d){0) =ui*S, dt{u*d){0) ~U2*S. Hence by uniqueness, u^u*d. Summing up, 

d f f d d 

— £;c= dtupdx+ udtpdx = -—EF-—EK. 

□ 



4. Weak solutions 



Based on Theorem 3.1 we now prove existence of global weak solutions to the 
relativistic Vlasov-Klein-Gordon system. The following auxiliary result will allow 
us to derive a-priori bounds from conservation of energy, at least for appropriately 
restricted initial data: 
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Lemma 4.1. Let p g]3/2,oo] and l/p+l/q=l. In addition to the assumptions of 
Lemma 5. 1 we assume further that d>0 with Jd~l. Let {f,u) be a solution as 
obtained in Lemma S.i.. Then 



u{t,x) p{t,x)dx 



<C{f)\\d,uit)hEKity^^ t>0, 



where 

cCf) := 
Proof. For any i? > 0, 



1/2 



3-7/6 



(7 + 3 







/ fdv + 




l\v\<B. 


J\v\>R V / 



(g+3)/6 



1/9 



\hr''''\\ht- 



If we choose R such that the right hand side becomes minimal we obtain the estimate 

3/(9+3) 



p{t,x)<C,\\f{t,x,.)\\f'i+'^ (^J ,/TTWfdvJ 



where 



4^x1/(9+3)^3^3X9/(9+3) 

3 j 3 



We take this estimate to the power {q + 3)/{q + 2), integrate in x, apply Holder's 

o 

inequality, and observe that ||/(t)||p = ||/||p to obtain 

l|pWII(,+3)/(,+2)<C,||/||^/(^+3)i?KW^/('+^^ 

By assumption on p we have {q + 3)/{q + 2) > 6/5, and by interpolation, 

(4.1) 

By Holder's inequality. Young's inequality, and Sobolev's inequality. 



pudx 



<\\u{t)\\e\\p{t)^d\\e/5<^\\d,uit)h\\pme/5 



< 



p*dudx 

1 ^(.+3)/6|j;||(3-9)/6||;||,/6||5^~(,)||^^^(^)l/ 

'J? 



with S3 3(7r/2)^/'^, cf. |jT^, 8.3]; recall from the proof of Lemma 3.2 that u = u*d. 

□ 

We now proceed to the main result of the present paper. In its formulation we 
employ the space 



ri m6 

-^kin 



l/ll 



kin ■ 



^/l + \v^\f{x,v)\dxdv< 



00 



Theorem 4.2. Let f e LI-^^{R'^) n L'' (R^) for some pe[2,oo], f>0, uiei/i(M^) 
U2ii L'^{M.^), and assume that 



(9+3)/3 
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where l/p+l/q=l. Then there exists a global weak solution {f,u) of the relativis- 
tic Vlasov-Klein-Gordon system (1.1), (1-i), (1-5) with these initial data, more 
precisely, 

/eL°°([0,oo[xLP(M6)), ueL°°{[0,oo[xH\M.^j) 

with 

dtueL°°{[Q,oo[xL^{R^)) 
such that the following holds: 



(a) if,u) satisfies ^T^, ^llj), ( |77| ) in V'{]0,oo[xR^). 



(b) The mapping [0,oo[b t ^ {f{t),u{t),dtu{t)) £ L^{m^) x L^{1 

° o o 

weakly continuous with {f,u,dtu)(0) = {f,Ui,U2). 

(c) /(t)>0 a.e., ||/(t)||p< ll/llp, t>0, and dtp + dWj = in V'{]0,oo[xR^) 
where j{t,x) :— J vf{t,x,v)dv. The weak solution conserves mass: \\f{t)\\i — 
II /111 for a. a. t>0. 

Remark. Regardless of whether one considers data satisfying the restriction of the 
theorem as small or not it should be emphasized that this is not a small data result 
of the type known for example for the Vlasov-Maxwell system, cf. since it does 
not rely on the fields being small and corresponding dispersive effects of the free 
streaming Vlasov equation. In particular, it is worthwhile to note that there is no 
size restriction on the data for the Klein-Gordon field. 

I 1 ° ° n 

Proof of Theorem 4.2. By our assumption on /, C(/) < 2 where the left hand 
side is defined as in Lemma |0| . We choose es]0,l[ such that C(/)^ < 2£. Next we 
choose sequences (/„) inC;?°(R6), /„>0, (mi,„), (u2,„) inC~(K3) such that /„^/ 
in Liijj(R<')nLP(MS), ui,„->ui in H^{R^), uz'n^Ma in L'^{R^), and we require that 



supC(/„ 

nGN 



<2e. 



(4.2) 



Let dn e 
and define 5- 



be non-negative, even, of unit integral, and with supp(i„ C i?i/„(0), 
dn*dn- Denote by (/n,Mn) the regularized solution according to 
Theorem 3J_ with S„ replacing d and initial conditions 

o 

Let p„ J fndv and denote by u„ the solution to the Klein-Gordon equation with 
right hand side — ^ dn and initial data 

W„(0) =tii,„*d„, dtUn{0)=U2,n*dn. 

As in the proof of Lemma 3.2 it follows that u„ = it„ * (i„ . Using Lemma 4.1 and 
Cauchy's inequality we find that 



1 



1 



E > ||/«(i)||kin + 7Tl|5t"«W||2 + :^l|i«n(t)||2 + 7Tl|5xU«(t)| 



'C{h)\\d,U,M\2\\fn{m'^' 



2' 

1/; 

kin 



> (l-£)||/„(i)|| 



kin " 



-^\\dtUn{t)\ 



-^\\Un{t)\ 



1 1 supC(A)^ 

Z 4£ fcgpj 



Hence by conservation of energy and (4.2) there exists C>0 such that for all 
ngNandt>0, ||/„(t)||kin, \\dtu„{t)\\2, h„(t)|Ui<a Also, ||/„(i)l|2 = ||/n||2 is 
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bounded by interpolation. Since the bounds are independent of n and i, by ex- 
tracting appropriate subsequences (again denoted by the same indices) we obtain 
/ei°°([0,oo[;L2(R6))^ wGL°°([0,cx)[;i/i(R3)) with e i°°([0,cx)[;i2(R3)) ^^^^ 

/„ - /in L2([0,r]xR6)^ 
M„ u in L^{[Q,T]xM?), 
d^Un^d^u in i2([o,T]x]R3)^ 
dtUn — dtu in L^{[Q,T]y,W^) 
for aU r>0. Moreover, /e L°°([0,oo[;Lj^i„(R6)) since 



V 1 + v'^fdxdv — hm lim / \/l + v^ fndxdv <C . 



Since ||/n(i)||p= ll/nllp is bounded as well, (4.1) implies boundedness of ||pn(i)||6/5 



Hence without loss of generality we may assume that pn — p in L6/5([0,r]xR3) for 
any T. Moreover, ||un(0ll6 < C'||9a;Un(i)ll2, so we may finally suppose that Un^u 
in L^([0,T] xK'^) for any T. We claim that u satisfies the Klein- Gordon equation 
in the weak sense. Indeed, let (/3eC^(]0,oo[xR^). Then 

= y (^dfUn~ AUn + Un+ Pn* Sn) fdxdt 

{undt(p~UnA(p + Un(p+{pn*Sn)(p) dxdt 



{^udfip — uAip + uip + pip) dxdt 

which yields the claim; note that J{pn*Sn)'p — J Pn (<5n *p) and (5„ * (/s <y9 in ( 

Turning now to the Vlasov equation, we adapt an argument from p^ . Sec. 3]. 
For £>0 let Ce&C^i^), < < 1, = 1 on [e,T], supp^^ C [e/2,2T]. Then 

dt{fn£,e)+vdx{fn^s)^<iiVy{fn^eda:Un)+fnC 

Hence by the velocity-averaging lemma of Golse, Lions, Perthame and Sentis ((|; D) 
we have 

yR>Oy^peC^iBii{0)) 3C^C{R,i^) VneN,e>0: 
\jU-)U.,..v)ij{v)dv\&H^'S^^^^) 

||/Ce(0/n(.,.,«)^(«)d«ll/fi/4<C(||^,/„||i-H||/„65.zi„||i + ||/„eili)'^'. 
In the case p — oo the boundedness properties derived above already assure the 
boundednes of the sequence (/Cs( •)/«(• > ■ i'f^)V'(^^)'^^)jieN in i?^/'*(M x M'^). The case 
p < (X) is more delicate and will be treated below. Since the restriction operator from 
H^/''(IRxM^) to L^([0,r] xBfl(O)) is compact, by a diagonal sequence argument, 
for each ^l; £C^{R^) we may extract a subsequence (again denoted by (/n)„GN) 
independent of e := l/m, R' = m (to £ N) such that 



U{.,.,v)^Piv)dv^ f{.,.,v)^Jj{v)dv (4.3) 



inL^{]0,T[xBR,{0)). 



For p<oo we show the validity of (4.3) as follows: For 77>0 we define /9,,(t) 



T/{l + r]T). Then Pn°fn, in addition to satisfying Vlasov's equation (with force 
dxUn) and the same bounds as /„, is bounded in L°° , and we may use the previous 
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case. In the limit 77^ we obtain L^-convergence in (4.3) which together with the 
uniform integrabihty of {{J f^''{.,.,v)tp{v)dv)'^} imphes the desired L^-convergence. 

We use this property to prove that / is a weak sohition of the Vlasov equation. 
For ^1 € (M) , (^2 , <^3 e (M3 ) , 

= y" JJ fn{t,X,v)(^ip[{t)Lp2{x)(p3{v)+V-dx(p2{x)ipi{t)ip3{v) 

- dxUnit,x) ■ dy(p3{v)(pi{t)ip2{x)) dvdxdt 



Choose a subsequence as above for 'ilj — dv'f2 and R' >0 such that supp(^2 Q Bri{0). 
Then 

J fn{., ■,v)dy(p3{v)dv^ J f{., .v)dy(p3{v)dv 

{[0,T]xBr,{ 

nally, 

f{t,X,v)(^ip[{t)ip2ix)(p3{v)+V-d^ip2ix)(piit)ip3{v) 



in LmO,T]xBR,{0)) and a,w„¥'i<^2 ^ 5,w(piv?2 m L^{{0,T)x Br,{0)). Thus, fi- 



= 



— dxu{t,x) ■ dyip3{v)ipi{t)(p2{x)j dvdxdt 

which estabhshes (a). 

As to (b), we integrate the Vlasov equation once with respect to t and define 

f{t):^f- v-dxf{s)ds- dxu{s)-dyf{s)ds 
Jo Jo 

(to be understood as an identity between distributions in x and v). Then, using the 
construction of / it is not hard to see that f{t) — f{t) in T>'{MP) for a.e. tE [0,oo[, 
that is, we can replace / by /. Moreover, /(O) = / and a straightforward application 
of the dominated convergence theorem shows that t'—^f{t) is weakly continuous. A 
similar argument applies to u. 

As to (c), we only prove that mass is conserved, the other assertions being 
standard. Let e>0 and choose R>0 such that J^^^^j^p{0,x)dx <e. Without 
loss of generality we can assume the same for pn{0), but since for the regular- 
ized problem the particles move along characteristics with speeds |w| < 1 we have 
J^^^yj^^j,pn{t,x)dx < e on any time interval [0,T]. Hence, 



p{t,x)dx > / p{t,x)dx— lim / pn{t,x)dx 

J\x\<R+T J\x\<R+T 

> lim / pn{t,x)dx — e=\\f\\i—e. 

n — >oo J 

o 

Since on the other hand J p{t, x)dx< \\f\\i and e is arbitrary this proves conservation 
of mass. □ 
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